DISTRIBUTION  STATEMENT  A.  Approved  for  public  release;  distribution  is  unlimited. 


Ocean  Acoustic  Propagation  and  Coherence  Studies 

Timothy  F.  Duda 

Applied  Ocean  Physics  and  Engineering  Department,  MS  1 1 
Woods  Hole  Oceanographic  Institution,  Woods  Hole,  MA  02543 
phone:  (508)  289-2495  fax:  (508)  457-2194  email:  tduda@whoi.edu 

Award  Number:  N00014-14-1-0223 
http://www.whoi.edu/people/tduda 


LONG-TERM  GOALS 

The  long-term  objective  is  to  develop  knowledge  of  the  factors  controlling  the  mean  and  the  variability 
of  acoustic  fields  in  the  shallow  ocean  environment  in  the  acoustic  frequency  band  50  Hz  to  3000  Hz. 
Propagation  variability  is  an  inescapable  complicating  factor  for  both  active  and  passive  sonar  systems, 
and  for  underwater  acoustic  communications.  Reliable  predictions  of  temporal  and  spatial  variability 
of  received  underwater  sound  can  improve  processing  and  handling  of  signals  of  interest,  for  example 
the  remediation  of  signal  degradation,  the  exploitation  of  intermittent  sonic  information,  or  the 
exploitation  of  favorable  propagation  conditions. 

OBJECTIVES 

The  main  goal  is  to  discover,  analyze  and  describe  mean  properties  of  acoustic  fields  and  time -varying 
acoustic  effects  caused  by  prominent  deterministic  propagation-altering  events  and  phenomena. 
Example  acoustic  effects  would  be  variation  of  horizontal  coherence  length,  variation  of  signal  level, 
variation  of  array  gain  (coherent  spatial  field  summation),  and  properties  of  the  (space-  and  time- 
dependent)  spatial  covariance  matrices  of  acoustic  fields.  I  plan  to  study  specific  acoustic  effects 
caused  by  water-borne  features  that  the  community  believes  are  the  most  influential:  water-mass 
fronts,  internal-wave  packets,  beam-like  internal  tides,  and  highly  structured  seafloors.  These  features 
can  impart  non-Gaussian  form  to  the  acoustic  field  variations  that  arise  from  ever-present  stochastic 
features  and  may  be  Gaussian  in  the  absence  of  the  features  but  the  presence  of  only  Gaussian  sound- 
speed  variability.  Seafloor  features  alone  can  also  impart  spatially  non-Gaussian  spatial  field  variation, 
with  example  features  being  sand  ridges,  slopes,  valleys  and  canyons. 

A  secondary  goal  is  to  more  fully  describe  the  propagation-altering  oceanic  phenomena  themselves. 
Working  towards  this  goal  is  closely  coupled  with  acoustics  research,  because  we  first  need  to  know 
which  features  are  the  more  relevant  for  acoustics,  and  to  what  detail  we  need  to  know  them.  This 
environmental  research  would  be  undertaken  with  acoustics  as  the  main  motivating  factor. 

APPROACH 

Although  much  progress  has  been  made  in  describing  ocean  acoustic  variability,  numerous  factors 
motivate  continued  research,  among  them  the  high  complexity  of  the  ocean  environment,  the 
associated  complexity  of  the  propagated  sound,  and  the  lack  of  a  simple  conceptual  framework  for  the 


1 


basic  properties  of  shallow-water  propagation  regimes  that  depart  from  the  most  simple  range- 
independent  adiabatie  mode  situation.  The  researeh  direetion  is  toward  sueh  a  framework,  to  be 
exploited  in  the  use  of  underwater  sound  in  shallow  water.  Our  knowledge  of  aeoustie  field  patterns  in 
shallow  water,  building  bloek  towards  a  useful  shallow-water  aeoustie  framework,  has  made  strong 
advanees  during  the  last  20  years  through  field  studies  and  theoretieal  studies.  We  wish  to  eontinue  this 
progression,  establishing  the  relative  real-world  importanee  of  various  fiuetuation-relevant  oeean 
features,  how  the  features  interaet,  and  how  the  sound  behaves  in  the  interaeting  features.  In  this 
approaeh,  models  of  the  environment  are  developed,  through  whieh  sound  propagation  is  simulated 
using  both  theoretieal  and  purely  eomputational  methods,  ineluding  our  3D  parabolie  equations  solvers 
(Lin  et  ah,  2013).  The  methodology  ineludes  building  environmental  models  of  hierarehieal 
eomplexity,  from  whieh  the  salient  aeoustie  proeesses  and  meehanisms  ean  be  identified  and 
quantified. 

The  aeoustie  simulation  tools  and  the  environmental  models  are  both  key  to  the  studies.  The  models  of 
the  variable  waveguide  features  ean  be  both  idealized  and  realistie,  with  idealized  models  used  to  study 
parameter  dependence  and  uneover  basie  physies,  and  realistic  models  used  for  comparing  computed 
and/or  theoretical  results  with  experiment.  Lying  between  these  two  activities  is  an  implieit 
eomparison  of  results  from  the  idealized  and  realistie  models,  whieh  may  answer  the  question  ean 
insights  obtained  from  the  simple  models  be  used  to  understand  results  that  arise  from  the  realistie 
models? 

An  insight  that  we  have  aequired  this  year  is  that  it  may  be  benefieial  to  examine  the  raw  output  of  the 
models,  for  instanee  phasor  behavior  in  time  and  spaee.  Software  that  ean  objectively  expose  and 
categorize  propagation  behavior  may  also  use  be  useful  to  eonstruct,  but  simply  looking  direetly  at 
transient  aeoustie  effects  can  be  valuable,  for  instanee  as  apparent  in  Badiey  et  al.  (2005),  Duda  et  al. 
(2012)  and  Badiey  et  al.  (2014)  where  strong  deterministic  acoustic  behaviors  that  are  non-Gaussian 
(at  least  over  the  studied  time  intervals)  are  examined  and  explained. 

The  oeean  models  and  the  aeoustie  simulation  tools  will  be  run  by  the  PI  and  by  Arthur  Newhall,  in 
eollaboration  with  other  members  of  the  WHOI  Oeean  Aeousties  and  Signals  Laboratory.  Five 
aeoustie  studies  are  planned; 


1 .  Canyon  and  slope  aeousties:  Identify  purely  geometrieally  eontrolled  (no  oeean  dynamics) 
vs.  time-dependent  aeoustie  field  features  and  statisties 

2.  Rough-seafloor  shelf  aeousties:  Determine  whether  findings  eoneerning  aeoustie  field 
eomplexity  from  water-eolumn  features  are  robust  to  bathymetric  and  sub-bottom 
uncertainty 

3.  Compare  the  roles  of  bathymetry  fiuetuations,  sub-bottom  properties  variations,  in  the 
presenee  of  internal  waves 

4.  Comparison  of  simulated  aeoustie  fields  resulting  from  realistie  oeean  environment  inputs 
and  those  resulting  from  idealized  oeean  environment  inputs 

5.  Evaluate  parameters  assoeiated  with  array  signal  proeessing  and  adaptive  signal  proeessing 

6.  Canyon  or  shallow-water  experiment  planning;  (2015-16  for  an  antieipated  2017  or  beyond 
experiment) 
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WORK  COMPLETED 


In  Study  1 ,  hierarchies  of  simulations  are  planned  for  a  few  bathymetric  geometries,  starting  with  a 
uniform  sound  speed  in  water  above  relief,  adding  refraction  from  background  sound-speed  profdes, 
adding  spatially-dependent  stratification,  and  finally  temporally  and  spatially-dependent  stratification. 
We  have  begun  with  the  most  complex.  The  4D-varying  ocean  fields  are  taken  from  the  S.  Calif  Bight 
internal-tide  response  simulations  (Ponte  and  Cornuelle,  2013).  Simplifying  the  environment  as 
planned  to  analyze  complexity  will  be  done  sequentially. 

Simulations  of  acoustic  fields  at  300  Hz  in  Coronado  Canyon  (Normark  et  ah,  2009)  have  been  added 
to  the  700-Hz  fields  examined  last  year.  The  change  gives  broader  area  coverage  (modeling  capability 
scales  with  acoustic  wavelength).  Figure  1  shows  bathymetry  and  output.  Simulations  are  at  one-hour 
intervals.  The  ocean  structure  contains  mainly  tidal  frequency  fluctuations,  which  are  resolves  at  this 
sampling  rate.  Some  flow  energy  may  be  present  at  higher  harmonics  of  the  tidal  frequencies,  as  well 
as  residual  subtidal  flows. 


E 


X  (km) 

HR03  250  m 


X  (km) 


-15 

-20 

-25 

-30 


-15 

-20 

-25 

-30 


Figure  1:  Snapshots  of  sound  intensity  (transmission  loss,  dB  scale)  at  two  depths  for  3D  parabolic 
equation-derived  simulation  outputs  for  300-Hz  sound  going  down  Coronado  Canyon  from  left  to 
right  are  shown.  The  source  is  40  m  deep  at  fx,y]=f0,0]).  White  lines  are  depth  contours  with  100  m 

increment.  The  500-m  contour  is  thicker. 
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Statistics  of  the  hourly  300-Hz  outputs  have  been  examined.  As  in  our  prior  work,  horizontal 
variability  in  the  acoustie  field  'F(x,y,z)  was  quantified  in  terms  of  the  eorrelation  function 


Rn{x) 


(fi^*(xo)fi^(xo  +  x)) 


where  the  averaging  (braekets)  ean  be  made  by  time  averaging  for  stationary  ergodie  processes,  and 
may  inelude  spatially  averaging  overxo  for  fields  with  homogeneous  statisties  (i.e.  Rn  is  not  a  funetion 
of  xo,  just  as  it  as  written).  Rn  was  examined  in  detail  by  Duda  et  al.  (2012)  for  field  data  from  a  plaee 
where  internal  wave  effeets  were  strong  and  the  seabed  was  shallow  and  relatively  flat.  Note  that 
f?v(0)=l.  The  point  x  =  L  where  Rn{L)  =  exp(-l)  is  referred  to  as  the  horizontal  correlation  length  scale 
(or  eoherenee  length),  L.  The  lagged  produet  quantifies  field  variability,  usually  dominated  by  phase 
differenees  for  typical  ocean  aeoustie  fields.  The  produet  behavior  is  linked  to  array  gain  for  additive 
beamformer  processors,  whieh  use  sums  of  spatially  lagged  signals  rather  than  the  products.  Gain  and 
correlation  can  be  expressed  in  terms  of  one  another  for  certain  funetional  forms  of  Rn  (e.g.  Cox, 
1973).  We  have  not  proven  it,  but  there  is  not  likely  to  be  a  general  relationship  between  the  additive 
and  multiplieative  behaviors  of  fields  with  uneonstrained  Rn  strueture. 

The  field  'F  and  the  eomputation  deserve  some  discussion.  To  analyze,  we  first  beamform  the  signals, 
as  was  done  with  the  SW06  data  (Duda  et  al.  2012),  eliminating  phase  wrap.  Figure  2  shows 
beamformed  output  (100-m  arrays)  at  two  loeations  in  the  eanyon  area.  This  is  an  example  of  the  near¬ 
raw  output  that  we  find  useful  to  examine. 
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Figure  2:  Beamformed  (array-steered)  output  for  two  array  locations  at  60-m  depth  in  the  canyon 
are  shown  (level  in  dB).  On  the  left  an  arrival  from  one  direction,  pulsating  in  time,  is  seen  for  a 
location  along  the  centerline  of  the  x,y  calculation  grid  over  the  canyon,.  On  the  right,  output  from 
a  site  close  to  the  canyon  edge  shows  multiple  arrivals,  one  at  high  angle  that  pulsates  and  two  or 
more  interfering  at  small  angle  (directed  at  the  source). 
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For  a  fluctuating  process,  the  time  mean  is  often  subtraeted  from  eaeh  point,  and  the  residual  forms  the 
field  to  be  analyzed.  These  fields  ean  be  analyzed  to  form  the  spatially-lagged  eovarianee  matrix  (part 
of  the  spatiotemporal  eorrelation  matrix  R  for  horizontal  apertures,  Johnson  and  Dudgeon,  1993).  This 
matrix  yields  for  example,  the  prineipal  components  (empirical  orthogonal  eigenfunctions)  that 
efficiently  describe  the  fluetuations.  This  would  provide  the  spatial  coherenee  behavior  of  the 
fiuetuating  field,  but  the  mean  values  may  need  to  be  retained  to  make  the  best  use  of  an  extended 
array  for  deteetion  or  elassifieation  problems.  We  are  working  to  examine  relationships  between 
principal  component  behavior,  eorrelation  behavior,  and  array  gain  behavior. 

This  suggests  a  seeond  option  where  field-measured  mean  values,  or  mean  values  compared  to  an  a 
priori  model  are  used.  A  simple  a  priori  model  is  the  eurved  or  spherieal  wavefront  used  for  nearfield 
beamforming.  A  recent  paper  from  our  group  studied  apparent  aeoustie  field  deeorrelation  from  oeean 
features  sueh  as  front  intervening  between  a  sound  souree  and  an  array  loeation,  whieh  will  impart  an 
unknown  but  steady  strueture  to  'F  (Lyneh  et  ah,  2014).  If  known,  the  feature  and  the  structure  could 
be  aecounted  for  to  aehieve  gain  over  a  very  large  array  aperture  (i.e.  creating  a  very  large  effective  L), 
but  currently  most  ocean  features  are  not  well  enough  mapped  to  be  aeeounted  for,  and  aetual  L  and 
effeetive  L  are  both  finite. 

We  eontinue  to  look  into  the  implieations  of  partial  field  eoherenee.  The  array  gain  AG  for  an  N- 
element  line  array,  the  power  of  the  beamformed  signal  divided  by  the  power  of  the  beamformed  noise, 
divided  the  signal  to  noise  ratio  for  a  single  sensor,  is  10  logic  {N^ /N)  for  a  perfectly  coherent  signal 
and  incoherent  noise.  This  ean  be  written  as  AG  =  SNRarray/SNRsensor,  where  SNRsensor  = 
Ssensor/Nsensor  for  a  single  sensor.  For  partial-array  coherenee,  the  coherence  length  can  give  AG  for 
an  arbitrarily  long  array  by  substituting  N  +(L/h,yf-  {L/^y)  for  in  the  numerator,  where  Ay  is  the 
inter-element  spaeing  in  meters,  so  that  AG=10  logic  {^^[{L/AyY -  (T/Ay)]/A0-  The  quantity  L/Ay  is  the 
number  of  effeetively  coherent  sensors  in  the  array.  Noise  is  eonsidered  uncorrelated  along  the  array  in 
this  model.  This  estimate  of  AG  is  useful  because  it  may  be  able  to  distinquish  between  non-monotonie 
Rn  behavior  of  multipath  interferenee,  and  monotonically  deereasing  eoherenee,  by  comparison  with 
direct  AG  computation.  In  log  units,  SNRarray  =  AG  +  Ssensor  -  Nsensor.  SNRarray  is  tied  elosely  to 
the  deteetability  of  a  partieular  signal  of  interest.  As  long  as  AG  +  Ssensor  exeeeds  Nsensor,  the 
ambient  noise  level,  the  system  is  usable.  Plots  of  results  from  this  caleulation  are  presented  in  the 
Results  section. 

RESULTS 

Figure  3  plots  correlation  and  array  gain  analysis  for  the  300-Hz  canyon  example  of  Figure  1.  The 
correlation  length  scales  at  60-m  depth  for  signals  steered  toward  the  source  with  100-m  long  synthetie 
arrays  (34  elements)  are  shown.  Also  shown  is  signal  exeess  Ssensor-Nsensor  at  all  positions,  again  for 
60-m  depth,  with  source  level  150  dB  and  Nsensor  95  dB.  Lastly,  a  panel  shows  SNRarray  with  AG  as 
in  the  model,  from  the  depicted  L. 

The  array  gain  can  of  course  be  eomputed  direetly  from  the  simulation  with  standard  delay  and  sum 
beamforming,  to  check  whether  the  model  that  uses  the  behavior  of  the  eross-produet  (the  L  estimate) 
to  estimate  AG  is  aecurate.  Agreement  between  the  two  would  indieate  that  the  Rn  is  generally 
monotonically  decreasing  and  that  the  scale  L  gives  usable  array  size.  Figure  4  shows  a  good 
eomparison  along  the  line  x  =  7700  m,  for  24  hours.  A  tidal  frequeney  oscillation  of  AG  is  evident  in 
the  figure.  The  output  of  the  beamforming  is  variable  over  time,  but  there  are  persistent  patterns. 
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Figure  3:  3D  parabolic  equation-derived  simulation  outputs  for  700-Hz  sound  going  down 
Coronado  Canyon  from  left  to  right  (150  dB  source  is  40  m  deep  at  the  [x,y]  origin.  Black  lines  are 
depth  contours  as  in  Figure  1.  The  top  shows  the  horizontal  correlation  length  L  for  lag  intervals  in 
they  direction,  calculated  at  60  m  depth,  computed  for  beams  directed  at  the  source.  The  calculation 
is  made  with  a  100-m  aperture,  so  computed  L  cannot  exceed  100  m.  The  center  plot  shows  the 
average  signal  excess  for  sensors  on  the  imaginary  arrays  with  respect  to  a  noise  level  of  95  dB.  The 
lowest  plot  shows  signal  excess  with  respect  to  the  same  noise  level  after  including  estimated  array 

gain  as  discussed  in  the  text. 


Figure  5  shows  beamformer  snapshots  for  two  lines  of  array  loeations,  again  at  60-m  depth,  along  with 
the  24-hr  means  of  the  outputs.  The  snapshots  depart  in  detail  from  the  mean  (not  yet  quantified)  but 
eapture  the  major  patterns  of  the  signal  direetionality  and  strength. 


The  24-hr  time  series  along  the  34-element  arrays  have  been  deeomposed  using  prineipal  eomponent 
analysis  (empirieal  modes).  For  the  domain  shown  in  Figure  3,  the  statistics  of  energy  explained  by 
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mode  two  (the  most  energetie)  divided  by  the  energy  explained  by  mode  one  (most  energetie)  is 
quite  variable.  At  10%  of  the  sites  the  ratio  is  greater  than  0.27.  At  2%  the  ratio  is  greater  than  .52.  At 
1%  it  is  greater  than  0.59.  A  single  dominant  mode  would  indicate  that  multipath  arrivals  fluctuate  in 
unison;  this  might  be  expected  for  tidal  effects  on  sound  speed  that  are  spatially  coherent  over  the 
modeled  domain.  Departures  from  this  are  under  examination.  We  are  now  looking  into  the  correlation 
of  significant  multimode  fluctuation  and  array  gain  performance.  It  is  possible  that  they  are 
independent.  Preliminary  analysis  of  related  properties,  direct-beam  power  versus  mode-one 
fluctuation  content,  is  shown  in  Figure  6.  Single  eof-mode  fluctuation  behavior  is  loosely  tied  to  more 
power  in  the  main  arrival  coming  straight  from  the  source  at  this  range,  (x  of  7700  m). 


Figure  4:  The  L-based  array  gain  estimate  (left)  which  results  from  analysis  of  beamformed  data 
snapshots,  is  compared  to  direct  estimates  of  array  gain  made  by  summing  the  signals  along  arrays, 
after  beamforming  (right).  Both  calculations  treat  noise  as  uncorrelated  between  sensors.  The 
computation  is  made  hourly  along  a  line  in  the  y-direction  atx=  7.7  km.  Gain  is  higher  over  deep 
water,  with  the  bathymetry  plotted  in  Figure  3.  Gain  with  full  coherence  of  34  elements  is  15  dB 


Figure  5:  (left)  24-hour  mean  beamformer  output  is  shown  for  a  lien  of  arrays  placed  atx  =  4200  m 
(see  Figures  1  and  3).  Zero  degrees  point  at  the  source.  So-called  out-of-plane  sound  at  higher 
angles  is  seen  in  the  average,  and  is  thus  time-coherent.  (2"‘'from  left)  One  snapshot  of  the 
beamformer  output  is  shown.  The  peak  levels  exceed  the  mean  levels.  (2'“' from  right)  The  24-hour 
mean  from  the  line  at  x  =  7700  m  is  shown.  The  peak  level  is  reduced  at  this  range  when  compared 
to  the  closer  range,  and  many  locations  show  very  little  coherent  arriving  sound,  (right)  A  snapshot 
at  7700  m  range  shows  some  transient  arrivals  that  are  not  evident  in  the  mean. 


1 


Direct  beam  vs  eof  fraction,  x=7700  m 


Figure  6:  Direct-beam  power  (at  the  angle  pointing  from  array  to  source)  versus  mode-one  principal 

component  fractional  fluctuation  energy,  is  shown. 

IMPACT/APPLICATIONS 

The  ultimate  performanee  limitations  of  marine  aeoustic  systems  are  set  by  the  properties  of  the  sound 
signals  that  propagate  within  the  fluctuating  three-dimensional  oceanic  environment.  Knowledge  of  the 
features  that  impact  sound  fields  and  how  the  features  affect  the  sound  can  guide  design  and  use  of 
acoustic  systems. 

RELATED  PROJECTS 

The  PI  also  has  a  MURI  grant  related  managed  by  ONR.  The  project  is  working  toward  coupled 
modeling  of  the  ocean  environment  and  acoustics  behavior  within  the  environment.  Information  can  be 
found  at  http://www.whoi.edu/sites/IODA.  The  PI  also  has  a  National  Science  Foundation  physical 
Oceanography  grant,  with  Co-PIs  Y.-.T.  Lin  and  Pierre  Lermusiaux,  to  study  ocean  internal  tide 
dynamics. 
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